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Abstract With an aim to examine whether the predicted solar gravitational 
redshift can be observationally confirmed under the influence of the convective 
Doppler shift due to granular motions, we attempted measuring the absolute 
spectral line-shifts on a large number of points over the solar disk based on an 
extensive set of 5188-5212 A region spectra taken through an iodine-cell with 
the Solar Domeless Telescope at Hida Observatory. The resulting heliocentric 
line shifts at the meridian line (where no rotational shift exists), which were 
derived by finding the best-fit parameterized model spectrum with the observed 
spectrum and corrected for the earth's motion, turned out to be weakly position- 
dependent as R::!-|-400 m near the disk center and increasing toward the 
limb up to sa-|-600 m (both with a standard deviation of cr w 100 m s""'^). 
■ Interestingly, this trend tends to disappear when the convective shift due to 

granular motions (w —300 m s^^ at the disk center and increasing toward the 
limb; simulated based on the two-component model along with the empirical 
center-to-limb variation) is subtracted, finally resulting in the averaged shift of 
698 m s^^ (<T = 113 m s~^). Considering the ambiguities involved in the absolute 
wavelength calibration or in the correction due to convective Doppler shifts (at 
least several tens m s~^, or more likely up to < 100 m s~^), we may regard that 
this value is well consistent with the expected gravitational redshift of 633 m s^^ . 

Keywords: Gravitational Redshift; Spectrum, Visible; Velocity Fields, Photo- 
sphere 



1. Introduction 

Ever since Einstein (1911) predicted based on his general theory of relativity 
that the wavelengths of solar spectrum lines should be shifted redward by a 
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proportion of 2 x 10^® 600 m s^^ in the velocity scale) as compared to the 
laboratory wavelengths, a number of solar physicists have tried to confirm this 
"solar gravitational redshift" observationally. Although very few people nowa- 
days doubt the practical validity of this prediction from general relativity, given 
the remarkable success in detecting this effect to a high precision from ground- 
based experiments, such as the well-known Pound-Rebka experiment using the 
Mossbauer effect (Pound and Rebka, 1960) or the experiment using two hydrogen 
masers on the ground and on a spacecraft (Vessot et at, 1980), efforts of solar 
astronomers toward confirming this gravitational redshift in the spectrum of 
the Sun have been non-terminating and continuing up to now. Apart from the 
historical studies in the early time of first 50 years which are reviewed in Forbes 
(1961); several notable trials using modern techniques and/or observational data 
of higher quality have been done over the last 50 years, as briefly summarized 
in Table 1. 

While the authors of these studies in Table 1 reported that they somehow 
detected the predicted redshift even to a quantitatively satisfactory degree, their 
results are not yet necessarily convincing, since they do not seem to have properly 
taken into consideration the effect of convcctive shifts due to granular motions 
and the closely related problem of line-shift determinations (i.e., definition of 
centroid wavelength) in case of asymmetric line profiles. Besides, most of these 
investigations are based on observations only on the disk center (or several 
disk points at best). Given that the well-known center-to-limb variation of the 
spectral line shift (i.e., the shift increasing redward toward the limb, often called 
"solar limb effect," first discovered by Halm (1907)) is regarded to be related with 
the photospheric velocity fields, it is highly desirable to study absolute line-shifts 
on a large number of disk points while considering the convective shifts (varying 
from point to point) at the same time. Namely, if the resulting residual shifts 
obtained from many points are confirmed to be nearly position- independent, it 
may suggest that shifts due to granular motions are successfully removed. 

Recently, we carried out an extensive spectroscopic study of solar differential 
rotation based on the spectra collected at a large number (« 3000) of points 
over the solar disk while applying the iodine-cell technique for Doppler-shift 
measurements (Takeda and Ueno, 2011; hereinafter referred to as Paper I). Since 
these data are suitable for the purpose mentioned above, we decided to conduct 
absolute Doppler shift measurements by using these spectra, while paying special 
attention to their subset observed on the solar meridian where the rotational 
shift is irrelevant, in order to examine how and whether the solar gravitational 
redshift may be successfully confirmed by eliminating the shifts of velocity-field 
origin. This is the aim of this investigation. 

The remainder of this paper is organized as follows. The observational data 
are explained in Section 2, and Section 3 describes our adopted method of 
analysis for absolute line-shift measurements. The resulting solutions and the 
simulation for the convective shifts to be subtracted are described in Sections 4 
and 5, respectively, based on which the final results for the solar gravitational 
redshift are discussed in Section 6. The conclusion is summarized in Section 
7. Finally, a supplementary section is presented as an appendix, where the limb 
effect empirically studied by Balthasar (1984) is reformulated for our application. 
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2. Observational Data 

We use the same observational data as adopted in Paper I, which were obtained 
on 20 and 21 July 2010 by using the 60 cm Domelcss Solar Telescope (DST) with 
the horizontal spectrograph of the Hida Observatory of Kyoto University (Nakai 
and Hattori, 1985) and the iodine cell placed in front of the spectrograph slit. 
One sequence is made up of 13 scanning observations with a step of Ar = R/12 
from the limb to the disk center at a fixed direction angle {6), and this sequence 
was repeated at 48 ^-values (with an increment of 7.°5), resulting in 13x48 
points on the disk as a total. The whole set was similarly conducted on each 
of the two days but with different slit setting; i.e., E-W direction (20 July) 
and N~S direction (21 July). Since three spectra (at different positions along the 
direction of the slit) were extracted from one observation, 3744 (= 13 x 48 x 3 x 2) 
spectra in the 5188-5212 A region (with reference I2 molecular lines imprinted) 
covering the entire solar disk were obtained. See Section 2 of Paper I for more 
details on the observational data. Prior to the analysis, a tentative wavelength 
calibration was done by using the representative solar lines, and an appropriate 
A (wavelengtlQ) vs. x (pixel number) relation was assigned to each spectrum 
frame. 

We examined the instrumental profile (-Pa) by fitting the observed lamp+l2 
spectrum {D'^^) with the modeled spectrum (-D^"" = Tx (E) P\) where T\ is 
the transmission function of I2 vapor and <Si means convolution. Regarding Tx , 
we adopt the spectrum (S/N 1000 with the resolving power of i? ss 400000) 
recorded for the I2 cell at Lick Observatory in the 5000-6300 A region by using 
the Fourier Transform Spectrometer at the McMath telescope of the National 
Solar Observatories at Kitt Peak (Marcy and Butler, 1992), which was kindly 
provided by Dr. G. W. Marcy (see also Takeda et ai, 2002). 

The instrumental profile was modeled by the sum of nine multiple Gaussian 
functions as P{v) (x '^'l^_^o,iexp[{v — i)"^] where v (in km s~^) is measured 
from the line center. Assuming oq = 1, we determined the 8 coefficients (a_4, 
a_3, a-2, fl-i; fli, 0,2, 03, 04) by requiring the best fit between the modeled and 
observed spectra in the 5202-5208 A region (Figure la). For a reference purpose, 
we also tried fitting with the Gaussian instrumental profile (cx exp[—{v/vni)'^]) by 
finding the best-match solution of the parameter Vm- The shape of the resulting 
P{v) is depicted in Figure lb, where we can see that the instrumental profile 
corresponding to our spectra (somewhat more rounded than the Gaussian) does 
not show any significant asymmetry that may affect absolute determinations of 
line shifts. The value of FWHM (« 5 km s~^) also suggests that the wavelength 
resolution is i? w 60000. 



^Unless otherwise explicitly stated, we mean "air wavelength in the standard condition" by 
this term. 
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3. Method of Absolute Line-Shift Determination 
3.1. Spectrum Modeling 

The iodine-cell technique has been extensively applied for precise Doppler shift 
determinations so far {e.g., spectroscopic search for exoplanets) since the pio- 
neering study of Marcy and Butler (1992) who showed how to accomplish high 
precision Doppler analysis based on a stellar spectrum with many superimposed 
I2 lines. However, its application has been limited essentially to "relative" mea- 
surements (i.e., with respect to a template spectrum whose absolute wavelength 
scale does not need to be precisely calibrated) so far. Since the determination of 
"absolute" shift is under question in this study, we have to modify Marcy and 
Butler's (1992) procedure (see Section 4 of their paper) appropriately. Namely, 
theoretically simulated spectrum with firmly established absolute wavelength 
scale must be used as the stellar template spectrum to be referenced. 

Generally, in analogy with the profile modeling conventionally used in stellar 
physics, the theoretical solar intensity spectrum with imprinted I2 lines at a disk 
point of direction cosine /i[= -^/l — (x^ + y'^)/R'^] may be expressed as follows: 



/th(A; fi) = fc{[^intr(A; fi; ^, Ai, A2, A3, . . .) ® Af^t(A; C)] • Ti. (A)} ® P(A), (1) 



where k is an arbitrary normalization factor, /i,itr(A;/x) is the intrinsic intensity 
at fi computed from a solar atmospheric model for given ^ (microturbulence) and 
Ai, A2, A3, . . . (elemental abundances), Afnit(A; fi; C) is the broadening function 
corresponding to macroturbulence Cj ^i^ (A) is the reference I2 spectrum men- 
tioned in Section 2, P{X) is the instrumental profile (Figure lb), and "(X)" denotes 
convolution (while "•" is usual multiplication). Note that such a line-broadening 
formulation based on the dichotomous micro-/macro-turbulence is nothing but 
a too simple modeling of actual velocity fields. Given that information for an 
adequate functional form of Af,„t is lacking, we move Mmt(A) to the outside of the 
square bracket in Equation (1) and further approximate Mmt(A; fi; C)®P(A) by a 
single Gaussian macrobroadening function G{X; wm), which is parameterized by 
the macrobroadening parameter vm (e-folding width) as G{v) oc exp(— u^/w^). 
Thus, we adopt the following model for line-spectrum simulation. 



where AAs(= A^pp — A^''^) and AAi2(= AjPp — X\f) are the apparent wavelength 
shifts of solar lines and I2 lines relative to the reference wavelength system {i.e., 
wavelengths given Table 2 for the former or in Table 3 for the latter), respectively. 

The parameters characterizing Ith (fc, AAg, AA12 , um, ^i, ^2, ^3, • ■ •) ^-'"^ 
termined by requiring that x^(= /[/th(A*: A) — /obs(M, A)]^dA is minimized (c/. 
Takeda (1995) for detailed descriptions of the numerical procedure). When the 
solutions of AAg and AAij have been established, the absolute line shift is 
obtained as 



/th(/^, A; k, AAs, AXi^,VM,Ai,A2, A3, . . .) 

= fc{[/intr(A + AAs; m; Ai, A2, A3, . . .) • Ti,{X + AAjJ} ® G(A; vm). 



(2) 



AAabs = AA, - AAl 



(3) 
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since AA12 is the error (due to instrument instability, etc.) to be subtracted from 

3.2. Atomic Line Data for Spectrum Synthesis 

As mentioned in Section 3.1, spectrum synthesis calculation for simulating the 
intrinsic specific intensity spectrum (/mtr) is important, which needs to be well 
reproduce the observed solar spectrum by adequately adjusting the elemental 
abundances and the macrobroadening parameter. Regarding the calculation of 
^intr(A,/i), we use Kurucz's (1993) WIDTH9 program, which was modified by 
Y. Takeda to enable spectrum synthesis by including many lines, along with 
Kurucz's (1993) ATLAS9 solar photospheric model (LTE, plane-parallel model). 
In addition, since the role of microturbulence is comparatively less significant 
in line broadening, ^ = 1 km is assumed as fixed. As to the atomic line 
data (wavelengths, oscillator strengths, etc.) necessary for this simulation, we 
adopt the data compiled by Kurucz and Bell (1995). Since trial test calculations 
suggested the necessity of some adjustments, the following modifications were ap- 
plied to in the \oggf values of four lines. Fe i 5197.929 (logg/ = -1.64 -1.34; 
increased by 0.3 dex), Fe i 5203.691 {loggf = -3.093 -9.99; neglected), Ni i 
5203.788 (loggf = -1.163 -> -9.99; neglected), and Fe i 5209.884 [loggf = 
—3.26 — > —3.46; decreased by 0.2 dex). Besides, since the weak depression at 
A w 5202A could not be reproduced (i.e., relevant line is missing in Kurucz and 
Bell's database), we neglected (masked) this local feature in evaluation. 

The finally adopted data for the important lines of appreciable strengths 
are summarized in Table 2, and how the observed disk-center spectrum (taken 
from the FTS solar spectrum atlas; Neckcl 1994, 1999jl can be fitted with the 
theoretical calculated /i = 1 intensity spectrum by using these atomic data (with 
appropriately adjusted elemental abundances and the broadening parameter) is 
displayed in Figure 2, where the lines given in Table 2 are also identified. 

3.3. Strengths and Formation Depths of Spectral Lines 

It is necessary to have knowledge on the strength and related quantities {e.g., 
line-forming depth) of spectral lines existing in the relevant wavelength range 
(5188-5212 A), since the convective line shift (its correction is important as 
will be discussed in Section 5) appreciably depends on such individual line 
properties. For this purpose, by using the solutions of the elemental abundances 
derived as by-products of the synthetic spectrum fitting in Section 3.2 (Figure 
2), we inversely computed EW (equivalent width at the disk center) for all lines 
presented in Table 2, while assuming that each line is single and isolated (i.e., 
irrespective of whether they arc actually blended or not). Besides, we computed 
two kinds of line-formation depth; one is the mean formation depth averaged 
over the profile ((logr)) and the other is the formation depth at the line-center 



^This spectrum atlas is available at ( ftp.hs.tini-hambtirg.de/pub/outgoing/FTS-atlas ). 
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(log To), which are defined as follows. 

n \ /-Ra logT5ooo(TA = l)<iA 

" 7i?IdA ' 

where R\{= 1 — I\/Icont) is the line depth of the emergent intensity profile, 
T5ooo('''A = 1) is the continuum optical depth at 5000 A corresponding to ta = 1 
at wavelength A within the profile, and the integration range is over the whole 
profile. 

log To = logr5ooo(TAo = 1) (5) 

where r^o is the optical depth at the line center. (All these quantities are com- 
puted for the /K = 1 ray at the disk center.) The resulting values of EW, (logr), 
and log To are given in Table 2. 

3.4. Wavelength Scale of I2 Reference Spectrum 

Regarding the reference spectrum of I2 molecular lines (T\), we adopt Marcy 
and Butler (1992)'s very high-resolution transmittance spectrum of Lick I2 cell 
already mentioned Section 2. This standard I2 spectrum in the relevant wave- 
length region, which we also used to fit the spectrum of our I2 cell (Figure la), 
is given in Table 3 (electronic data) . Though the original data are represented in 
terms of the vacuum wavelength scale (Avac), we converted it into air wavelength 
(Aair; which is usually used in optical spectroscopy such as in Kurucz and Bell's 
atomic-line database we adopted) as 

Aair Avac/^air ■ (6) 

Here, riair is the refractive index of air in the standard condition represented as 
a function of the wavenumber k cm^^ 10^/Avac(A)], for which we adopted 
the formula given in Cox (2000): 

nair = l.+6432.8xl0"^+2949810./(146.xl0^-A:2)+25540./(41.xl0^-A:2). (7) 

While this Lick I2 spectrum has so far been successfully used for very precisely 
detecting "relative" variations of radial velocity intended for exoplanet search 
or stellar seismology, it is not clear whether its "absolute" wavelength scale is 
sufficiently reliable. Since this point is of primary importance in this study, we 
examined its precision by comparing it with the theoretical I2 spectrum by using 
"lodineSpec" program (Knockel et ai, 2004), which can model the wavelengths 
of the rovibronic structure of the I2 B-X spectrum with very high precision (its 
uncertainty is on the order of only w 30 MHz at « 5200 A, which means a 
precision of « 5 x 10~® or « 15 m s~^; cf. their Figure 3), though the intensity 
distribution is not so precisely reproduced as wavelengths. The comparisons are 
displayed in Figures 3a (wide view) and 3b (magnified view of the 5196-5197 A 
region), where the vacuum wavelength scale (Avac) is used in the abscissa. We 
notice from Figure 3b that the Lick I2 spectrum appears to be slightly shifted 
redwards relative to the theoretical one. This is quantitatively confirmed by 
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a cross-correlation analysis (Figure 3c), showing that the average shift in this 
5196-5197 A bin amounts to 0.0028 A or w 160 m s^^. We similarly examined 
all twenty-four 1 A-bins in the 5188-5212 A region and the histogram of the 
resulting shifts are depicted in Figure 3d, from which we conclude that the 
absolute wavelength scale of the Lick I2 spectrum is systematically redder and 
should be corrected blueward by —167 m s~^ (simple average of 24 values; with 
a standard deviation of 42 m s~^). 

Besides, given that we transformed the vacuum wavelengths of reference I2 
spectrum into the air wavelengths by applying Equations (6) and (7), the air 
wavelengths of the atomic lines we adopted for simulating the solar spectrum 
have to be essentially on the same system. For the purpose of checking the 
consistency, we examined the Avac/'^air ratios (or the refractive index in the air; 
ukb) for all 719 lines (5188-5212 A region) included in Kurucz and BcU (1995) 
as 

^KB = (-^uppcr £^lower)/(10«Aair), (8) 

where i?uppcr and -Eiowcr arc the excitation potentials of lower and upper levels 
in cm~^, respectively, and Aair is the air wavelength in A (all these three kinds 
of data are presented in Kurucz and Bell's database). The comparisons of such 
evaluated tike's and the corresponding formula values (^formula) resulting from 
Equation (7) are depicted in Figure 4. We can recognize from this figure that 
the consistency is satisfactorily good and the average of the relative difference 
[(("■KB — "■formula)/'^)] IS Only —4.47 X 10~^ (with a standard deviation of cr = 
6.46X 10-8) or -13 m s"! (cr 19 m a-^)E Thus, even though small, we formally 
apply a correction of —(—13) m s^^ to the air wavelengths of I2 lines (c/. given 
3rd column of Table 3) in order to maintain consistency with the wavelengths 
of atomic lines given in Table 2. (Regarding the sign of the correction, since 

;^formula ^ ;^KB becauSC of HkB < ^formula, WC havC tO shift A^'?™™''' slightly 

redward to make it consistent with Kurucz and Bell's system.) 

Consequently, we conclude that the A}^^ values given in Table 2, which we 
adopt as the reference (air) wavelengths of I2 lines, should be corrected blueward 
by (5A}^^ = —154 m s^^ (net combination of —167 m s"^ and -1-13 m s"^). More 
precisely, since the wavelengths of I2 lines directly affect AJ®^ of AAij (= Aj^*^'" — 
Xlf), it is included in AAabs as the term of positive sign (-l-Aj^^; cf. Equation 
(3)), which means that AAabs should be corrected by +SX\f ~ —154 m s~-^. 
In practice, the whole fitting analysis is conducted by literally using the (air) 
wavelengths of I2 spectrum given Tables 3, and this correction is applied to the 
results at the last stage. 



^This dispersion (a) of !^ ±20 m s^^ is reasonably attributed to the fact that the wavelengths 
of Kurucz and Bell's (1995) line data are presented only to the third decimal of A (cf. Table 
2), which means that the intrinsic uncertainty amounting up to < 0.0005 A or < 30 m is 
inevitable in our analysis. 
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4. Results 



Although our primary interest is only on the meridian points {x w 0) unaffected 
by the rotational velocity, we carried out our absolute wavelength-shift analysis 
for all the available data points over the entire solar disk, except for those at the 
limb because the solutions there were found to be comparatively uncertain and 
unstable {i.e., the data at ri2 were not used; cf. Figure 2 in Paper I). 

Our spectrum covering 5188-5212 A was divided into four segments of 5- 
6 A [regions (a), (b), (c), and (d), cf. Table 2], and AA^^j^g was determined for 
each segment i (i = 1, 2, 3, 4; each corresponding to a, b, c, d, respectively) by 
applying the procedure described in Section 3.1. An example of how the observed 
solar-|-l2 spectrum matches the modeled spectrum corresponding to the best-fit 
solutions of the parameters is demonstrated in Figure 5. Then, the (raw) radial 
velocity of the observed point on the solar surface (V^^^^^) and its probable error 
(e) were computed as 

V, = cAAlbJA, - 154 (m s'^), (9) 



yr.^ EE (10) 



and 

~N{N - 1) 



V (v- - yI■aw^2 



where we applied the correction of —154 m s^^ due to the absolute wavelength 
scale problem of the adopted I2 spectrum (Section 3.3), c is the speed of light, 
and N = 4. Further, we added the heliocentric correction for the earth's motion 
j-^hei. ^j- ggction 4.1 in Paper I) and obtained V^°^ in the heliocentric system as 



Figure 6 demonstrates, for the meridian case, how the results of V[ obtained at 
different observational times become mutually consistent by applying the A'^''' 
correction. 

The distribution of the probable error (e) for each measurement is shown in 
Figure 7, which is directly comparable with Figure 9 of Paper I. As can be seen 
from Figure 7a, e mostly ranges between sa 50 m s~^ and 100 m s^^; that is, by 
w 2-3 times larger than the case of relative measurements in Paper I, reflecting 
the difficulty of absolute measurements. We do not sec any clear dependence 
of e upon the position (r) on the disk, except for a slight decrease near to the 
limb (Figure 7b). The results of our analysis [observational time, {r,6) as well 
as {x,y) coordinates, VJ."^™, A^'^', and V^^^] are presented in Table 4 (electronic 
data) for the meridian case of our concern, which will be discussed in Section 6. 
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5. Simulation of Convective Blue-Shift 

5.1. Line Shift due to Convective Granular Motions 

Although we have derived the absolute line-shift reduced to the rest 

system of the Sun at each point on the disk, it alone is still insufhcient for 
our main purpose of quantitatively confirming the gravitational redshift, since 
the systematic shift due to convective velocity fields has still to be subtracted. 

A hot up-rising bright bubble (granular cell) originating from the deeper con- 
vective zone gradually decelerates its speed as it moves upward while dissipating 
its energy, and eventually stops and turns into cool downflows (dark intergranular 
lane). Since the contribution of the brighter former is larger than the darker latter 
in the integrated spectrum, the net effect generally makes a blue shift as seen 
from the observer. Since the extent of velocity contrast is depth-dependent (i.e., 
becoming less significant toward higher layers), such a blue shift tends to be 
more appreciable for weak deep-forming lines than strong high-forming lines, as 
seen in Figure 2 of AUendc Pricto and Garcia Lopez (1998). For the same reason, 
the off-center region (forming deeper) of a line profile is more blue-shifted than 
the region of line-center core (forming higher) , which results in an asymmetric 
profile with a "C" -shaped bisector (see, e.g., Dravins, 1982). 

To make things further more complicated, not only the vertical velocity field 
but also the horizontal field has to be considered when studying the variation 
of the convective shifts over the disk. This results in a fiatter decreasing rate 
of blue shift (from disk center toward the limb) compared to the case of the 
vertical field only; actually, even a peak of blue shift is occasionally observed not 
the disk center but at off-center around fi w 0.8 (Schrotcr, 1959; Beckers and 
Nelson, 1978; Balthasar, 1984, 1985). 

Considering this complexity, we decided to adopt the following strategy: As 
far as the convective shift at the disk center (/x = 1) is concerned, we can evaluate 
it for any relevant line to a sufficient precision, thanks to the availability of a 
well-tuned two-component model. Meanwhile, since this kind of line shift at 
positions off the center (0 < ^ < 1) is too difficult to simulate, we invoke the 
empirical work of Balthasar (1984), who extensively investigated the relative 
/x-dependence of the shifts for many spectral lines of different properties and 
classified them in terms of line- formation depth, which we reformulated in a more 
easy-to-use manner (c/. the Appendix). Accordingly, as successively described 
in the following sections, we (1) first simulate the absolute line-shifts at the 
disk center (/.i = 1) for representative lines (Section 5.2), and (2) then combine 
them with the relative center-to-limb trends empirically elucidated by Balthasar 
(1984) (Section 5.3). 

5.2. Absolute Shift at the Disk Center 

For the purpose of simulating the disk-center intensity profile under convective 
velocity fields, we adopt the two-component model proposed by Borrero and 
Bellot Rubio (2002). Their model consists of two kinds of plane- parallel atmo- 
spheric models (with only vertical velocity fields) corresponding to hot/bright 
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granular (g) and cool/dark intergranular (i) region. The details are presented 
in Appendix A in their paper, where T (temperature), Po (electron pressure), 
Pg (gas pressure), p (density), and v (vertical velocity) arc given as functions of 
depth (along with the depth-independent microturbulcncc and macroturbulence) 
for each of the two models. The structures of T(t5ooo) and w(t5ooo) for each model 
are depicted in Figure 8a. 

By using the WIDTH9 program (as in Section 3) along with these models, 
we computed the emergent specific intensity profile /^(A) and /'(A) at the disk 
center, from which the final spectrum is eventually obtained as the combination 
of these two: 

^comb(^) = (1 - a)/g(A) + aPiX), (13) 

where a is the filling factor assumed to be 0.24 according to Figure 3 of Borrero 
and Bellot Rubio (2002). 

In order to check our computational procedure, we simulated the /i = 1 
profiles of 22 Fe i lines given in Table 1 of Borrero and Bellot Rubio (2002), and 
confirmed that mosl|f| of the resulting line-shifts satisfactorily reproduced those 
of the observed line profiles of Neckel's (1994, 1999) solar disk-center spectrum 
atlas. Figure 8b shows an example of such simulated I^, P, and (each being 

normalized with respect to the continuum level of its own) for the representative 
Fe I 5194.941 line. 

Such computed P°'^^\X) is further convolved with the instrumental profile 
(P; cf. Section 2), resulting in an asymmetric broadened profile corresponding 
to our instrumental condition. Profile examples before and after LP. convolution 
are displayed in Figure 8c again for the case of Fe i 5194.942. 

The final step is to determine the line shift from such simulated (asymmetric) 
profile by applying the same procedure described in Section 3; i.e., fitting with 
the classically modeled symmetric profile while adjusting the line shift {5X or 
5V) along with the broadening width and the elemental abundance. This process 
is demonstrated in Figure 8d, where we can see that the match is satisfactory. 

We applied the whole above-mentioned procedure to 18 lines of appreciable 
strengths, which we selected by the criterion of (logr) < — 1. Such determined 
line shifts (SVi) at the disk center (/z = 1) for these lines are given in the last 
column of Table 2, where we can see that SVi is in the range of —310 m s~^ < 
SVi < -230 m s-^l 



More precisely, the predicted shifts of the 3 lines (Fe I 6574.2285, Fe I 6581.2100, and 
Fe I 6739.5220) out of 22 lines appear to be only slightly shifted bluewards compared to 
the observation; otherwise, the profile/ wavelength match for the remaining lines was found to 
be very good. 

^ It should be stressed that we adopted the fitting- based line-shift (SVat), while the line shift 
((5Vniin) defined in terms of the wavelength of minimum intensity has been more commonly 
used so far. We also determined i5Vniin by interpolating the lowest three profile points by a 
quadratic function (see, e.g., AUende Prieto and Garcia Lopez, 1998) for all these 18 lines 
in the present case, and found that the difference SV^nin — SVat tends to be EW-depedent, 
which varies from -60 m s"! (EW 70 mA) to -|-120 m s"! (EW 230 mA) and can be 
approximated by the linear-regression relation <5Vmin - "SVflt (ms"!) = -129 + 1.02EW (mA). 
Since the typical extent of this difference is several tens of m s~^, we may state that the results 
arc not seriously influenced by the definition of line shift. 
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5.3. Center-to-Limb Variation of 5V 

We are now ready to derive the ^-dependent convective shift, (5V^, by combining 
the disk-center SVi computed in Section 5.1 and the center-to-limb variation 
formulated in the Appendix based on the results of Balthasar (1984). 

For a line k, we first evaluate logTg (Balthasar's core-forming depth) from 
logTQ (our line-center-forming depth) by solving Equation (23). Then, the coef- 
ficients and c'^ in the velocity unit can be derived from Equations (21), which 
are sufficient with Equation (19) to establish the relative /^-dependence of SV as 

SV'^{^,) - SV'il) - A4(/.) - A4(1) ^b'^{l-pi) + c\l - (14) 

Therefore, combining Equation (14) with SVi which was already derived in 
Section 5.2, we can express the absolute convective shift for line k at any fj, as 

sv^- = SV,' + b'{i + c'ii - (i5) 

Such obtained position/angle-dcpcndcnt (5V^' results for the 18 lines arc shown 

in Figure 11a {5V vs. fjL) and lib {5V vs. •\/l — /i^), where we can see that the 
tendencies are nearly similar to each other, except for two conspicuously strong 
Cr I fines (at 5206.038 and 5208.419 A) of low excitation (xiow « 1 cV). 

As the final step, we average these (5V^ results over 18 lines in order to obtain 
the mean relation, while adopting 

w'' (X 1/{t') = 10-<'°s^'> (16) 

as the weigh10 for line k. 

Accordingly, wc finally obtain 

{5V^) - {6V^) + {b){l-^Ji) + {c){l-^Jif 

= -278.5 + 79.9(1 -/i)+ 422.3(1 -/i)^ (m s"^), (17) 

where (• • •) is the w'^-weighted average over 18 lines. Such obtained {5V^), which 
is also depicted in Figures 11a and lib by the thick red line, will be used in 
Section 6 for correcting V^"^ for the convective shift. 



® This choice of weight (ex l/(r)) is for the following reason. The precision of radial velocity 
determination is higher not only for stronger (deeper) but also for sharper (narrower) lines, 
which means that the radial velocity is mainly determined by sharp saturated lines at the fiat 
part of the curve of growth. Meanwhile, "very" strong lines with appreciable damping wings 
do not necessarily contribute so significantly because of their less sharpness. We point out 
here that the extent of saturation (i.e., both of strength and sharpness), to which more weight 
should be given, is most well described by 1/(t) because it increases as a line gets stronger 
from the linear to the fiat part of the curve of growth, while it then turns to decrease as a line 
further becomes very stronger and comes into the damping part of the curve of growth, which 
is the characteristic just what we want. (In contrast, 1/ro is ever increasing monotonically 
with the growth of line strength.) Figure lid describes well this situation. 
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6. Discussion 

We now discuss the results of absolute line shifts obtained in Section 4, while 
taking into account the convcctive-shift correction estimated in Section 5. In 
the following description, x and y denote the Cartesian coordinates on the 
solar disk (c/. Figure 2a in Paper I). The relations between V^''^ obtained for 
all the observed points and each of the positional coordinates are displayed in 
Figures 10a [V^^"^ vs. y) and 10c {V^''^ vs. x), where the exhibited trends can be 
reasonably understood by the Doppler shift due to solar rotation {V^'^^ oc x; c/. 
Section 4.1 in Paper I). A rough inspection of these two figures suffices to realize 
that V^'^'^ distributes around w +500 m s~^ at cc « (meridian line) or near 
to the disk center (where the data points are densely confined) , which suggests 
that V^°^ approaches the expected gravitational shifl|j of +633 m s^^ where the 
rotational Doppler shift is negligible. 

In order to study this situation more in detail, we pay attention to the selected 
dataset corresponding to the observations on the meridian line (a; w 0), which 
are given in Table 4. The V^"^^ vs. y correlation for this subset is shown in Figure 
10b, where we can confirm that V^''^ values confine around ~ +500 m s~^ (with 
a standard deviation of a ~ 100 m s^^) as mentioned above. In addition, we can 
recognize a tendency of slight increase of 1/.*^°' with |y|; i.e., sa+400 m s~^ at 
the disk center and increasing toward the limb up to ~+600 m s~^. (Actually, 
such a change of V,.'^"' from the limb to the center on the meridian line is already 
observed in Figure 6.) This trend is nothing but the "limb effect" already men- 
tioned in (see, e.g., Figure 2 of Dravins, 1982). Interestingly, when we overplot 
the {5V) vs. y relation according to Equation (17) [describing {5V) as a function 
of jjL (= yjl — {y/RY] derived in Section 5, its tendency is found to be similar 
to the behavior of V^'^''' (c/. the solid line in Figure 10b). Thus, 1^'^°' is corrected 
for {5V) as 

V^^^-^°" = V^^^ - {5V), (18) 

and the resulting ]/hei,corr plotted against y in Figure lOd. We can see in this 
figure that the slight curvature observed in Figure 10b is more or less removed by 
this correction (especially in the northern hemisphere, though somewhat over- 
corrected in the southern hemisphere) and X/'i^i.corr ^^-^ nearly j/-independent 
along the meridian line, which is just what we intended to achieve {i.e., to 
consistently detect the gravitational redshift at a number of different points by 
eliminating the position-dependent limb effect; c/. Section 1). The distribution 
histogram of l/hoi,corr g^^^^ ^]^^ ^^hoi,corr^ (locallv averaged value) vs. y relation 
are shown in Figures lOe and lOf, respectively^ The average of all the "(Aei.corr 
(meridian) data turned out to be ((1/^'^°'^'^°'''^)) ~ 698 m s~^ (with a standard 
deviation of a = 113 m s^^). 



'^This value of 633 m s ^ corresponds to an observer on the earth while the redshift is 636 m s ^ 
for the observer at infinity. 

® It is interesting to note that our {V^'^^"'^"^'' ) of 696 m s~^ at the disk center (c/. Figure lOf) is 
comparable to 684 m s~^ derived by Balthasar (1985) as the asymptotic value for the absolute 
disk-center line-shift at the strong-line limit of logTB — > — oo (where the convective shift is 
expected to be insignificant). 
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Here, it may worthwhile to recall that several sources of errors are involved 
ill our results: Each individual absolute V[ measurement is accompanied by a 
probable error of sa 50-100 m s~^ (Figure 7). Also, the randomness of granular 
velocity fields (represented by "wiggly-line" character seen in high-resolution 
slit-spectrograms) produces an appreciable dispersion of Vi- (even though it must 
have been somewhat smoothed because spatial information over 42" has been 
averaged in our spectrum; cf. Section 2 in Paper I), which is estimated in the 
present case to be on the order of « 100 m s~^ (Figures lOe and lOf). Yet, since 
these are random errors, they could (in principle) be reduced in the averaged 
{{Vr}) by summing up as many data as possible; e.g.^ averaging of all 144 data on 
the meridian would result in a reduction by a factor of w 10 (sa -\/l44). Wc thus 
estimate that the error component in ((V,-)) due to such random sources would 
have been reduced to the « 10-20 m s~^ level. Actually, however, systematic 
errors must be more important, which are hard to evaluate. For example, while 
we corrected the absolute wavelength-scale of Lick I2 spectrum by 162 m s~^, 
this value itself already has an uncertainty with a magnitude of several tens 
m s~^ (Figure 3d). Besides, although our convective-shift corrections turned 
out seemingly successful, we have no idea about how much systematic errors 
arc involved in these, given that this correction differs from line to line with 
variations of ~ 50-100 m s~^ {cf. Figure 9). Considering all what described 
above, we would modestly state that our final {{V^'^'^^''^°")) value suffers inevitable 
ambiguities of no less than several tens m s~^ or maybe amounting up to < 
100 m s"^ 

Accordingly, while the resuhing ((V;'^'='''=°")) of 698 m s ^ is slightly larger 
than the predicted gravitational redshift of 633 m s^^ by -1-65 m s~^, this 
difference is not meaningful in view of the expected uncertainties of several 
tens m s~^ to < 100 m s~^ as estimated above. We thus conclude that the solar 
gravitational redshift has been quantitatively confirmed in a reasonable manner 
by our absolute line-shift analysis using the iodine-cell technique, in combination 
with the simulated convective shift corrections. 

7. Conclusion 

Motivated by the situation that various past studies of solar gravitational redshift 
are not convincing because few of them seem to have properly considered the 
Doppler shifts due to inhomogeneous convective motions, we decided to exam- 
ine whether the predicted solar gravitational redshift can be observationally 
confirmed by removing the appropriately simulated convective shift. 

For this purpose, we attempted measuring the absolute spectral line-shifts on 
a large number of points over the solar disk based on an extensive set of 5188- 
5212 A region spectra taken through an iodine-cell with the Solar Domeless 
Telescope at Hida Observatory, where we paid attention only to the selected 
dataset observed on the solar meridian where the rotational Doppler shift is 
irrelevant. 

In analogy with the analysis technique for precise "relative" shift measure- 
ments (from l2+object spectra) elaborated by Marcy and Butler (1992), we 
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developed a method for measuring "absolute" shifts by using the theoretical solar 
intensity spectrum template, which is computed from a model atmosphere and 
parameterized by elemental abundances along with the broadening parameter; 
and the absolute shift can be established by finding the best-fit solution of the 
parameterized model spectrum with the observed spectrum. 

The resulting heliocentric absolute shifts (V^.'^^') on the meridian line turned 
out to be weakly position-dependent as w+400 m near the disk center and 
increasing toward the limb up to w+600 m s^^ (both with a standard deviation 
of « 100 m s-i). 

We then simulated the convcctive blue shifts for representative 18 lines by us- 
ing Borrero and Bellot Rubio's (2002) two-component model for the disk center, 
combined them with the empirical relations of relative center-to-limb variation 
investigated by Balthasar (1984), and finally adopted the weighted-average over 
18 lines. The resulting shift is (SV) w —300 m s"-'^ at the disk center, gradually 
increases to « —100 m s~^ at \y/R\ ~ 0.9, and shows a suddenly increase up to 
« +200 m s"i at the limb {\y/R\ = 1). 

Interestingly, the position-dependence of V^°^ tends to be removed when the 
simulated convective blue-shift is subtracted, which finally makes the averaged 
shift for all the meridian data to be ((F,''^''™")) = 698 m s'^ (cr = 113 m s'^). 

Considering the ambiguities involved in the absolute wavelength calibration 
or in the correction due to convective Doppler shifts (at least several tens m s~^, 
or more likely up to < 100 m s~^), we may state that this value is well consistent 
with the expected value of -1-633 m s~^ (due to the gravitational potential on the 
solar surface) as predicted from the general relativity theory. This means that 
we could confirm the solar gravitational redshift in a quantitatively satisfactory 
level within the precision of our analysis. 

Acknowledgements Y. Takeda heartily thanks Dr. H. Knockel for kindly providing the 
"lodincSpcc" program, which turned out very helpful for checking the absolute wavelength 
scale of the reference I2 spectrum. 

Appendix 

Analytical Formulation of Balthasar's (1984) Empirical Limb-EfFect 

Balthasar (1984) carried out an empirical investigation on the center-to-limb 
variation of spectral-line shift on the solar disk and found that the trend of line- 
asymmetry as well as of line-shift depends on the formation heights of line cores 
(logre). He tried to fit the observed shifts of many spectral lines by a quadratic 
polynomial in terms of 1 — /i (/i = cos 6: direction cosine) as 

AAb (or Avb) = a + 6(1 - /i) + c(l - nf. (19) 

and determined the coefficients (a, b, c) for 6 different logre ranges (c/. Table IV 
therein), as shown in Figures lla~llc. For practical convenience, we derived the 
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analytical expression of a, b, and c with respect to logre by applying quadratic 
least-square fitting to these data as 

a (mA) = 0.254 + 0.6f9(logrB) +0.138(logTB)^ 

6(mA) = -7.03 + 0.00450(logrB) +0.8f9(logTB)^ (20) 

c(mA) = f8.4+1.90(logrB)-0.530(logrB)^ 

as depicted in these figures. Considering that Balthasar (1984) used 5000 A as 
the standard wavelength, we can convert these coefficients into velocity scale as 

a{ms-^) = 15.3 + 37.1(logrB) + 8.30(logTB)^ 

6(ms-i) = -421.8 + 2.70(logrB) +49.1(logrB)^ (21) 

c(ms"^) = 1100. + 113.9(logTB)-31.8(logrB)^ 

One problem is that Balthasar (1984) did not clarify the definition of his 
"line-core formation depth" (expressed as "Igr" therein but denoted here as 
logTB). We computed (logr) and logro (the mean formation depth and the line- 
center formation depth defined in Equations (4) and (5)) for all 143 lines given 
in Table I of Balthasar (1984) in the same way as we did in Section 3.3. These 
results are compared with his logrB values in Figures lid and lie, where we 
can see that his logrB is similar to our logro, though both are not in a linear 
relationship with each other. We found from least-squares analyses the following 
relations: 

(logr) = -0.0375 + 1.03(logrB) + 0.124(logrB)^ (22) 
logro = 0.0136+ 1.54(logrB) + 0.111(logrB)^ (23) 

In the practical application in Section 5, our logrg is converted into logrB by 
using Equation (23) for each line, from which the necessary coefficients (such as 
b or c) can be derived by applying Equation (21). 
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Table 1. Solar gravitational redshift determinations over the last 50 years (expressed by R = Aobs/^thc)- 
Authors: Blamont and Roddier (1961) 

Method; Sr I 4607 line with atomic-beam resonant scattering method 

Results: R 0.5 (si 5 mA/9.76 mA) [disk center], _R ss 1.0 (fa 10 mA/9.76 mA) [limb] 

Remark: Redshift of 2.4 mA due to Stark effect is subtracted 

Authors: Snider (1970) 

Method: K I 7699 line with atomic-beam resonant scattering method 
Results: R = 0.61 ± 0.06 (10 ± 1 mA/16.3 mA) [disk center] 

Authors: Snider (1972) 

Method: K I 7699 line with atomic-beam resonant scattering method 
Results: R = 1.01 ± 0.06 (16.4 ± 1 mA/16.3 mA) [disk center] 

Authors: Beckers (1977) 

Method: Measuring line-shifts of Zeeman-insensitive {g = 0) Ti I 5713 line in the umbrae of 7 sunspots 
Results: R = 0.97 ± 0.02 (613 ± 14 m s" 1/633 m s"!) 

Authors: LoPresto et al. (1980) 

Method: 738 Fe lines based on Pierce and Breckinridge's (1972) wavelength table (disk center) 

Results: R = 0.76 ± 0.24 (all 738 lines at 2930-5762 A), R = 0.97 ± 0.16 (only 74 strong lines at 3075-4958 A) 

Authors: LoPresto et al. (1991) 

Method: Wavelength shifts (at the profile core) of oxygen triplet lines at 7772-7775 A 

Results: R = 0.73 (460 m s-^/633 m s-^) [disk center], R = 0.99 (627 m 8-1/633 m s^^) [limb] 

Remark: Measurements on 7 different points suggest a tendency of increasing shifts toward the limb 

Authors: Cacciani et al. (2006) 

Method: Doppler shift measurements over a number of points over the disk by using a magneto-optical filter 

Results: R = 0.99 (625 m s^i/633 m s^^) [at the meridian line of zero-rotational velocity] 

Remark: Several assumptions are involved regarding other sources causing lino shifts (c/, section 6.1 therein) 

Note: The theoretical solar gravitational redshift of 633 m s" ^ (the value corresponding to 
observations from the earth) is adopted here. 
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Table 2. Atomic and relevant data of important lines at each region of spectrum fitting. 



Species 


(A) 
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(eV) 


log af 
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Gammas 


Gammaw 


EW 
(mA) 


(log 




log 


TO 


SVi 

(m s-1) 


Region 


(a): 5188.5-5193.6 A (Fc 


:, Ti, Ca, 


, Cr, 


Ni, and 


Nd abundances varied to 


adjust) 












Ti II 


5188.680 


1.582 


-1.210 




8.23 


-6.66 


-7.95 


91.7 


-1 


,52 


-2 


.42 


-249.6 


Ca I 


5188.844 


2.932 


-0.090 




8.42 


-4.18 


-7.31 


124.2 


-1 


.10 


-2 


.08 


-262.9 


Ti I 


5189.621 


2.239 


-1.017 




8.06 


-5.71 


-7.79 


3.7 


-0, 


,06 


-0, 


.07 




Fc I 


5191.455 


3.038 


-0.656 




8.00 


-5.47 


-7.62 


176.4 


-1 


,69 


-3, 


.68 


-294.5 


Cr I 


5192.002 


3.395 


-0.390 




7.91 


-5.32 


-7.62 


23.2 


-0, 


.33 


-0, 


.48 




Fc I 


5192.343 


2.998 


-0.521 




8.01 


-5.47 


-7.62 


204.3 


-1 


,64 


-3, 


.92 


-282.3 


Ni I 


5192.490 


3.699 


-1.512 




8.08 


-5.67 


-7.69 


26.7 


-0, 


.41 


-0, 


.60 




Nd II 


5192.614 


1.136 


0.290 




7.92 


-5.80 


-7.73 


8.7 


-0, 


,15 


-0, 


.22 




Ti I 


5192.969 


0.021 


-1.006 




6.78 


-6.32 


-7.84 


84.6 


-1 


,78 


-2 


.63 


-247.1 


Cr I 


5193.480 


3.422 


-0.720 




6.83 


-6.11 


-7.84 


11.9 


-0, 


,16 


-0, 


.23 




Region 


(b): 5193.6-1 


3199.0 A (Fe, Cr, Mn, Ni 


, Ti, and 


Y abundances varied to 
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Ti I 


5194.034 


2.103 


-0.560 




7.83 


-5.37 


-7.62 


10.2 


-0, 


.13 


-0, 


.19 




Fc I 


5194.941 


1.557 


-2.090 




6.29 


-6.20 


-7.87 


120.7 


-2, 


,26 


-3, 


.74 


-298.9 


Fc I 


5195.468 


4.220 


-0.002 




8.31 


-6.17 


-7.67 


102.9 


-1 


.43 


-2 


.40 


-277.0 


Fc I 


5196.077 


4.256 


-0.451 




8.33 


-5.16 


-7.67 


77.4 


-1 


,16 


-1 


.80 


-240.6 


Y II 


5196.422 


1.748 


-0.880 




7.92 


-5.93 


-7.77 


22.5 


-0, 


.37 


-0, 


.53 




Cr I 


5196.482 


3.449 


-0.270 




7.80 


-5.93 


-7.82 


21.1 


-0, 


,30 


-0, 


.44 




Cr I 


5196.590 


3.449 


-0.360 




7.83 


-5.94 


-7.82 


18.0 


-0, 


.25 


-0, 


.37 




Mn I 


5196.593 


3.135 


-0.930 




7.77 


-6.09 


-7.83 


9.5 


-0, 


.13 


-0, 


.19 




Ni I 


5197.157 


3.898 


-1.190 




8.14 


-5.59 


-7.81 


21.7 


-0, 


,32 


-0, 


.47 




Fe II 


5197.577 


3.230 


-2.100 




8.48 


-6.60 


-7.95 


86.8 


-1 


,17 


-1 


.93 


-250.3 


Fc I 


5197.929 


4.301 


-1.340 




8.33 


-4.77 


-7.77 


31.7 


-0, 


,47 


-0, 


.69 




Fc I 


5198.711 


2.223 


-2.135 




8.21 


-6.18 


-7.83 


87.5 


-1 


,67 


-2 


.55 


-253.3 


Region 


(c): 5199.0-5205.0 A (Cr, Fc, Y, 
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and Ti abundances 
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Fc I 


5199.531 


3.882 


-3.236 




8.01 


-5.47 


-7.69 


1.5 


-0, 


.03 


-0, 


.04 




Fc I 


5199.695 


3.546 


-3.248 




7.69 


-6.11 


-7.80 


2.9 


-0, 


,05 


-0, 


.06 




Nd II 


5200.121 


0.559 


-0.490 




7.92 


-5.89 


-7.75 


10.4 


-0, 


,18 


-0, 


.27 




Cr I 


5200.207 


3.385 


-0.660 




7.91 


-5.34 


-7.62 


14.7 


-0, 


.20 


-0, 


.29 




Y II 


5200.406 


0.992 


-0.570 




7.92 


-6.03 


-7.80 


31.5 


-0, 


.56 


-0, 


.83 




Ti I 
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Shown here arc the atomic data of the lines of appreciable strengths (sec also Figure 2), which were taken 
from Kurucz & Bell's (1995) compilation (with adjustments for 4 lines; cf. Section 3.2). Followed by first 
four self-explanatory columns, damping parameters are presented in columns 5—7: Gammar is the radiation 
damping width (s""*"), log-yj-g^^j. Gammas is the Stark damping width (s""*") per electron density (cm~^) at 
10'* K, log(7c/A'^G)- Gammaw is the van der Waals damping width (s""*") per hydrogen density (cm~^) at 10** 
K, log(7w /iVfj ) . In case where these damping parameters arc not given in Kurucz and Bell (1995), they were 
computed as default values in Kurucz's WIDTH program (Leushin and Topil'skaya, 1987). Columns 8—10 give 
the strength-related quantities relevant to the disk-center spectrum {cf. Section 3.3): the equivalent width 
inversely computed from the best-fit abundance solution {cf. Figure 2) under the assumption of isolated 
single line, the mean depth of line formation, and the line-center formation depth. Presented in column 11 
is the theoretically estimated convective blue-shift at ^ = 1 (only for the representative 18 lines; cf. Section 
5.2). 
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Figure 1. Determination of the instrumental profile {P\) by comparing the observed lamp+l2 
spectrum {D'^'') with the simulated spectrum (-D^™ = Tx Px, where is the transmission 
function of I2 vapor, (a) Upper spectra are the observed D'^" (open circles) with the best-fit 
Z)|'™ (solid line) overplotted, while the lower spectrum (solid line) is Tx- (b) Thick solid line is 
the finally obtained P{v) (the sum of nine off-center Gaussian functions; cf. Section 2), while 
the dashed line is the special case derived by assuming that P is the Gaussian form. The 
abscissa is in the velocity scale, and both profiles are normalized at the maximum peak to be 
unity. 
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Figure 2. Demonstration of how the solar intensity spectrum simulated at = 1 for each 
wavelength region (a— d) can be fitted with the observed solar disk-center spectrum (Neckcl, 
1994, 1999) by adjusting the abundances and the macrobroadoning parameter. Observed and 
theoretical spectrum are expressed by (red) open circles and (blue) solid lines, respectively. 
Identifications are given for the line features of appreciable strengths. 
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Figure 3. Comparison of the theoretical pure ^^^l2 absorption spectrum simulated by 
Knockel et ai.'s (2004) "lodineSpec" program with the observed I2 absorption spectrum (Lick 
I2 cell spectrum recorded with Kitt Peak FTS) adopted as the reference in this study: (a) wide 
view for the 5188-5212 A region and (b) magnified view for the 5196-5197 A region (note 
that vacuum wavelength scale is used in these two figures), (c) Result of cross-correlation 
analysis for the 5196-5197 A region shown in panel (b), where the abscissa corresponds to the 
wavelength correction to be applied to the observed spectrum to enable a match between the 
two. (d) Histogram showing the distribution of the corrections (expressed in the velocity scale) 
derived by cross-correlation analyses for each of the twenty-four 1 A bins in the 5188-5212 A 
region. 
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Figure 4. Wavelength-dependence of the relative difference between nformula a-nd tikb (two 
kinds of refractive indices) for all the spectral lines at 5188A < Aair < 5212A included in 
Kurucz and Bell's (1995) compilation, where nformula is the refractive index of air based on 
the standard formula given in Cox (2000) as given in Equation (7) while tikb is estimated 
from Kurucz and Bell's (1995) atomic data by Equation (8). 
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Figure 5. Graphical example [for the disk-center case at region (a)] of how the absolute 
wavelength shift can be determined by searching for the best-fit parameters such those ac- 
complishing the best-fit (simulated) I\ th [= G\ {Tx ■ Ix intr); c/. Equation (2)] with the 
observed Ix obs (observed solar intensity spectrum recorded through the I2 absorption cell. Top: 
comparison of Ix obs (red open circles) and Ix th (blue solid line). Middle: Tx (transmission 
spectrum of I2 vapor cell). Bottom: /A.intr (computed intrinsic solar intensity spectrum). 
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Figure 6. Runs of derived raw radial velocities (V/''™), heliocentric corrections (A*^"'), and 
heliocentric radial velocities (1/'^'^') with time, for four kinds of scanning observations along the 
solar meridian line relevant for this study, each of which were carried out from the limb to the 
disk center over 2 minutes. The results for the 2010 July 20 observations (east-west aligned 
slit) are shown in panels (a) (9 = 0°; northern meridian) and (b) (9 = 180°; southern meridian). 
Similarly, those for the 2010 July 21 observations (north-south aligned slit) are presented in 
panels (c) and (d). In each panel, the horizontal dashed line denotes the zero-level. 
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Figure 7. (a) Distribution histogram for the probable error (e) in Vr (where results all over 
the disk are combined), (b) Mean probable errors ((e); averaged at the same radius bin) plotted 
against the radius. 
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Figure 8. (a) r(r5ooo) (solid line) and ^(tsooo) (dashed line; positive for the downfiow) 
structures of two-component solar photospheric model (blue and red lines corresponding to 
granular and intergranular region, respectively) derived by Borrero and Bellot Rubio (2002), 
which we adopted for the line-profile simulation, (b) Illustration of how the realistic asymmetric 
profile (green) is produced as a combination of granular (blue) and intergranular (red) profiles 
(being weighted according to the filling factor and the continuum intensity), for the case of 
Fe I 5194.941 line at the disk-center, (c) Example of how the simulated profile (green solid 
line) in (b) is broadened to result in a realistic profile (black dashed line) corresponding to 
our observation by convolving the instrumental profile (c/. thick solid line in Figure lb), (d) 
Demonstration of how the asymmetric profile simulated in (c) (black-dashed line) can be 
fitted with the classically-modeled symmetric profile (pink solid line) by adjusting the profile 
parameters (Fe abundance, line broadening parameter um, and the line shift). The O — C 
residual is also shown at the bottom of the figure (with an offset of +0.1). 



SOLA: body.tex; 3 July 2012; 0:17; p. 25 



Y. Takcda and S. Ueno 




Figure 9. Angle/position-dependent convective line sliifts {SV) for 18 representative lines 
of appreciable strengths, which were simulated based on Borrero and Bellot Rubio's (2002) 
two-component model for the disk center (SVi) and combined with Balthasar's (1984) em- 
pirical SVfx — SVi relation (c/. Section 5.3 and Appendix), (a) 5V vs. cos9(= fi), (b) SV vs. 
sin6'(= y/R). The finally adopted relation (averaged over 18 lines) is depicted by the thick red 
line. 
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Figure 10. Behaviors of the V^°^ (absolute heliocentric Doppler shift in rn s^^) at various 
positions of the solar disk {x coordinate: E-W direction increasing towrard west, y coordinate: 
S— N direction increasing toward north); expressed in unit of the solar radius R). (a) V^'^"' vs. 
y for all measurements over the entire disk, (c) x vs. V^*""' for all measurements over the entire 
disk, (b) V^"^ vs. y only for the meridian measurements (x = 0). The positional-dependence 
of the simulated convective shift cf. Equation (17)], which is a function of y through 

the relation fj. = ^1 — (x^ -|- y'^^jB? = \/l — (s//-R)'^ is also depicted by the solid line, (d) 
y^hcljcorr^^ yliA _ ^^V^j)) VS. y Only for the meridian measurements {x = 0). (e) Histogram 
showing the distribution of (yhel^ where (V^*^"') is the average of data at 

almost the same meridian points (bars indicate the standard deviation). In panels (a)-(d), 
probable errors involved for each of the Vr measurements are indicated by thin ticks. The 
expected solar gravitational redshift at the earth (633 m s~^) is shown by the (red) arrow in 
all six panels. 
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Figure 11. Panels (a)-(c): Graphical representation of the depth-dependence of the coefR- 
cicnts (a, b, c) used for the polynomial fits of the limb effect, AA = a + 6(1 — fi) + c(l — ^j,)^, 
which were empirically determined by Balthasar (1984; cf. Table IV therein). The quadratic 
dependence of these coefficients in terms of logrg, which we derived from the least-squares 
fit are also depicted. Panels (d) and (e) show the relationships between the line-core-forming 
depth (logre) adopted by Balthasar (1984) and the two kinds of line-formation depths ((logr) 
and log to; i.e., mean depth of line formation and the line-center- forming depth) derived in 
this study. The least-squares fitted quadratic relations are also drawn. 
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